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CN ' We analyze the non-equilibrium dynamics of the O(N) $ 4 model in the 

large N limit with broken symmetry tree level potential and for states of large 
energy density. The dynamics is dramatically different when the energy den- 
sity is above the top of the tree level potential Vq than when it is below it. 
When the energy density is below Vq, we find that non-perturbative particle 
production through spinodal instabilities provides a dynamical mechanism 
for the Maxwell construction. The asymptotic values of the order parameter 
43 ■ only depend on the initial energy density and all values between the minima 

of the tree level potential are available, the asymptotic dynamical 'effective 
potential' is flat between the minima. When the energy density is larger than 
Vo, the evolution samples ergodically the broken symmetry states, as a conse- 
quence of non-perturbative particle production via parametric amplification. 
Furthermore, we examine the quantum dynamics of phase ordering into the 
broken symmetry phase and find novel scaling behavior of the correlation 
function. There is a crossover in the dynamical correlation length at a time 
scale t s ln(l/A). For t < t s the dynamical correlation length oc \ft 
and the evolution is dominated by linear instabilities and spinodal decom- 
position, whereas for t > t s the evolution is non-linear and dominated by 
the onset of non-equilibrium Bose-Einstein condensation of long-wavelength 
Goldstone bosons. In this regime a true scaling solution emerges with a non- 
perturbative anomalous scaling length dimension z = 1/2 and a dynamical 
correlation length oc (t — t s ). The equal time correlation function in this 
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scaling regime vanishes for r > 2(t — t s ) by causality. For t > t s phase or- 
dering proceeds by the formation of domains that grow at the speed of light, 
with non-perturbative condensates of Goldstone bosons and the equal time 
correlation function falls of as 1/r. A semiclassical but stochastic description 
emerges for time scales t > t s . Our results are compared to phase ordering in 
classical stochastic descriptions in condensed matter and cosmology. 

11.10.-z;11.15.Pg;11.30.Qc 
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I. INTRODUCTION 



Recent studies of non-equilibrium evolution of high energy density states have revealed 
novel and unexpected phenomena ranging from strong dissipative processes via particle 
production |],^| to novel aspects of symmetry breaking |||4j]. These non-equilibrium effects 
are crucial to understanding the end of inflation and the reheating problem @,[|,|| . They are 
also conjectured to be relevant for baryogenesis, the partial restoration of symmetries during 
phase transitions ]7|,|8| , for the formation of defects during non-equilibrium phase transitions 
HU and reheating dynamics in gauge and fermionic theories ||10|| . More recently a novel 



form of anomalous relaxation with non-universal critical exponents [IT] has been found in 
scalar field theories due to non-equilibrium quantum effects. 

One of the important aspects of these phenomena are their non-perturbative charac- 
ter. There are few consistent approaches that can be used for their description. While a 
classical approach has been advocated to address these issues [§,|7|-0| we have emphasized 
the application of the large N limit that captures both the quantum aspects as well as the 
emergence of classicality |p^,|6l,pT| . The large N limit in scalar field theories has a consis- 
tent implementation within the non-equilibrium quantum field theory formulation p|,p|JT^] 
and the renormalization aspects had been studied in detail fl3|jl4]] . Furthermore it provides 
a consistent, renormalizable framework to study non-perturbative phenomena that can be 
systematically improved [p],[2|.|T2l . 

In this article we focus on further new non-equilibrium aspects. First, we discuss how 
non-perturbative particle production due to spinodal instabilities gives rise to a dynamical 
version of the Maxwell construction of statistical mechanics. It is a well known result that 
the static effective potential becomes complex in the region of the phase diagram in which 
coexistence of phases can occur. The imaginary part has been deemed to arise from a 
failure of the loop expansion, but a more careful treatment reveals that this imaginary part 
is actually related to the decay rate of a quantum mechanically unstable state [jn|[nj. More 
precisely, the static effective potential approach is just not appropriate to study the dynamic 
of phase transitions. At the level of formal statistical mechanics, the mixed phase is described 
via a Maxwell construction that leads to an effective potential (free energy) that is flat in 
the coexistence region. The Maxwell constructed effective potential has been confirmed by 
a careful analysis using the 'exact' renormalization group |T7|^m[ of the equilibrium coarse 
grained free energy. 

Here we study the dynamical aspects of the Maxwell construction in the large N limit 
of a scalar field theory. We begin by reviewing the relevant features of the static effec- 
tive potential in the large N limit p0|-^1 to highlight the fact that the imaginary part 
of the effective potential is not an artifact of the loop expansion but a more fundamental 
shortcoming of a static and equilibrium description of the coexistence region in terms of 
an homogeneous order parameter. We then study in detail how non-equilibrium dynamics, 
and in particular, particle production due to spinodal instabilities, leads to a flat dynamic 
potential in the region of coexistence. We find that all expectation values of the scalar field 
between the minima of the tree level potential are available asymptotically and that the 
asymptotic value is reached with a vanishing effective mass for the expectation value. This 
is the dynamical equivalent of a 'flat dynamical potential' and the availability of all possible 
asymptotic expectation values between the minima is the dynamical equivalent of the lever 
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rule for coexistence of phases in the Maxwell construction. For a given initial expectation 
value, such that the total energy density, e, is smaller than the maximum of the tree level 
potential V(0), [here V($) is the potential for the 0(N) field $], we find that the asymptotic 
expectation value depends on the initial conditions only through the total energy as follows 



< oo) >= J —\m\ 



1 - 



V(0) 



where A and m are the renormalized quartic self-coupling and mass respectively and x a 
new dynamical anomalous exponent which is universal in the weak coupling limit and given 
by x = 0.25. 

In summary, the dynamical effective potential is always real and is flat for e < V(0). 
Two physically different situations appear here: 

a) when the energy density e is smaller than the maximum of the tree level potential, e < 
V(0) (in terms of renormalized mass and coupling), the non-equilibrium dynamical evolution 
leads to broken symmetry states and the Maxwell construction via non-perturbative particle 
production as a consequence of the spinodal instabilities. 

b) when e > V(0), in which case the energy density is larger than the maximum of 
the tree level potential, we find that the evolution is symmetric in that the time evolution 
of the expectation value of the scalar field samples the minima of the tree level potential 
symmetrically and ergodically with vanishing average over time scales much longer than 
the period. What we find is that this evolution is a consequence of the growth of the 
root mean square fluctuations due to parametric amplification which eventually induces a 
positive mass squared for the field. We show numerically that the expectation value relaxes 
to zero asymptotically transferring all of the energy to quantum fluctuations via parametric 
amplification in accord with the results in PJTT|. 

Our final topic concerns the quantum non-equilibrium dynamics of phase ordering in the 
O(N) model. This is an important topic and is the basis of the Kibble mechanism for defect 
formation as the early universe cools through phase transitions |53-E3J. Recently, there have 
been some attempts at constructing a non-equilibrium, real-time, description of the dynamics 
of this process |26|]. The early stages of the dynamics of phase ordering in quantum field 
theory models begin with the growth of spinodally unstable long-wavelength field modes and 
a dynamical correlation length emerges ||27|| . The dynamics of non-equilibrium fluctuations 
during the early spinodally unstable stage has been studied and clarified recently providing a 
deeper understanding of the quantum fluctuations that trigger the process of phase ordering 
p8| , |29[] . We go beyond the early and intermediate time regime which is dominated by 
the linear instabilities, into the highly non-linear regime which is dominated by non-linear 
I I] and that begins at the spinodal time t s ~ m _1 ln(l/A). This time scale 



resonances 



determines the emergence of a semiclassical but stochastic description which we discuss and 
clarify. After t s , semiclassical, large amplitude (oc l/VX) field configurations are represented 
in the quantum density matrix with unsupressed probability. 

We find that for t > t s a true scaling solution emerges for which the field acquires a non- 
perturbative anomalous scaling (length) dimension z — 1/2, a dynamical correlation length 
£(t) ~ (t — t s ) emerges. In and the equal-time correlation function vanishes by causality 
at distances r > 2(t — t s ). In the non-linear regime we find the onset of a novel form of 
Bose Einstein condensation in which the zero momentum mode of the quantum fluctuations 
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grows asymptotically linearly in time becoming macroscopically populated at long times. 
This condensation mechanism results in the equal time correlation function decreasing as 
1/r for r < 2(t — t s ) as a consequence of massless excitations. We find a crossover in the 
time dependence of the dynamical correlation length at the spinodal time t s : 

£(t) « Vi for t < t a ; f (t) &t-t s for t > t s . 

The process of phase ordering in this regime is described by the formation of domains 
that grow at the speed of light, inside of which there is a non-perturbative condensate of 
Goldstone bosons. 

In section II we review the construction of the effective potential in the large N limit 
and highlight several important features that are relevant for the dynamics namely 

• the regime of validity of the scalar field theory as an effective theory below the energy 
scale of the Landau pole ~ m e 167r / A , and the consistency of studying this theory for 
energy densities oc m 4 / A, 

• the complex result for the effective potential in the coexistence region and the failure 
of such effective potential to describe dynamical aspects. In the large N limit, this 
effective potential is exact and this complex value cannot be blamed on a failure of the 
loop expansion but is a consequence of more fundamental shortcomings of an effective 
potential description. 



In section III we set up our study of the non-equilibrium dynamics in terms of the quan- 
tum density matrix in the large N limit and its probabilistic interpretation. This formulation 
is the most enlightening to use use in order to discuss the emergence of classicality in terms 
of probability distributions. It also provides a natural framework in which to interpret the 
formation of semiclassical, large amplitude, field configurations that are represented in the 
quantum ensemble with unsuppressed probabilities. 

Section IV is devoted to a detailed study of the dynamics in a broken symmetry potential. 
Here we distinguish between symmetric and non-symmetric evolution and provide a simple 
criterion for ergodic symmetric evolution in a broken symmetry potential. We contrast the 
energy with the effective potential and show that particle production provides a dynamical 
mechanism for the Maxwell construction. Section VI deals with the non-equilibrium process 
of phase ordering both in the linear and in the non-linear regime. We then compare our 
results to those of phase ordering kinetics in condensed matter physics [p0| , pl| and the scaling 
solution found in radiation and matter dominated Friedmann-Robertson- Walker cosmologies 



for the non- linear sigma model p2| , p3| . We summarize our results and conclusions in section 
VII. 



II. THE EFFECTIVE POTENTIAL IN THE LARGE N LIMIT 

One of the main conclusions of our work is that the understanding of field theories to 
be found in the effective potential can be woefully lacking under certain circumstances. In 
order to better compare and contrast our dynamical results with what would have been 
expected from an effective potential analysis, and since there are some subtleties involved in 
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the discussion of the large N limit effective potential p0|- p^1 , we recapitulate this discussion 
below. 

There are two facets of the O(N) effective potential that we want to make contact with: 



• Scalar field theories can only be taken as effective field theories with a restricted regime 
of validity. The O(N) model is known to have a Landau ghost at an energy scale 
« m e 16n / A , with m ; A the renormalized mass and coupling respectively, and can 
only be used reliably below this energy scale. A fuller analysis of the O(N) effective 
potential in the large N limit determines the range of energy densities for which the 
theory can make sensible predictions. 

• The effective potential is complex in the spinodal region. This is what gives rise to the 
standard Maxwell construction. We want to understand the dynamics associated with 
this construction. 

We consider the O(N) vector model with Lagrangian density given by, 

C = ~ $(x, tf - ~ , t)f - V($(x, t)) 

^) = \rnt& + ^($r + ^ (2.1) 

To analyze the large N limit it proves convenient to introduce a Lagrange multiplier field 
a(x) [12|] and write the above Lagrangian density in the form 



£ = ~ !(£, tf - ~ [V$(x, t)} 2 + ||- (a(x, t) - 2 mjj) a{x, t)- 1 - tf a{x , t) 

The effective potential (V e ff) can be obtained in the N = oo limit in the following manner 
22|: introducing the expectation value of the scalar field as <p(x, 
the quadratic functional integral over the field <$>(x,t) can be performed exactly leading to 
an effective functional of the Lagrange multiplier a(x,t) and the expectation value <fr(x,tf. 
The leading order in the large N limit is obtained via the saddle point approximation in 
the functional variable a(x,t) and expansion around the saddle point generates a series in 
powers of 1/JV. Translational invariance dictates that the relevant saddle point solution 
must be a space-time constant. In the large N limit only coupling and mass renormalization 
are needed: 



Ao 
I 

Ao 




(2.2) 
(2.3) 



where A is an ultraviolet cutoff and the finite part in the renormalization of the coupling 
( p.3| ) has been conveniently chosen to simplify the expressions below. As a function of 
a(x) = M 2 , and <ft 2 , the effective potential at leading order in the large N limit is: [20 — 22| 
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eff\ 



M 



M 2 



M 2 - 2m 2 
A 



M 2 
32^ 



log 



e \m\ 



M 2 



(2.4) 



where m and A are the renormalized mass and coupling, and we have chosen \m\ 2 to be 
the renormalization scale. Extremizing with respect to the saddle point solution for the 



Lagrange multiplier, M we are led to the gap equation |20|-|22 

t2 



2(M 2 



m 



A 



M 2 , 
+ 16^ l0g 



\m\ 
~M 2 ~ 



(2.5) 



Shifting the Lagrange multiplier field a(x) by the saddle point value M 2 we identify M 2 
as the effective squared mass of the fluctuations of the <3> field. 

We see from eq. (|275|) that <p 2 (M 2 ) increases with M 2 for M 2 < \m\ 2 e 32n2 / x - 1 . 
The effective potential at the solution of the gap equation becomes 



V e// (0 2 (M 2 ),M 2 ) 



M 4 
16" 



I \m\ 2 \ 
A ' 4vr 2 \y/eM 2 ) 



+ 



log 



(2.6) 



It proves convenient to introduce the following dimensionless variables: 

M 2 



w 


= e -327r 2 /A 




32tt 2 


v(w) 


= 4 e 




m 4 




16tt 2 




= 1 12 e 




m z 



64 ^ 2 / A \/ e// (0 2 (M 2 ),M 2 ) 



in terms of which 



v(w) 



-w - + log w 



(2.7) 
(2.8) 

(2.9) 



In fig. Q we plot the dimensionless function v (w) vs. w. Furthermore at the solution of the 
gap equation we find the relation 



ip 2 (w) 



32vr 



sign(m 2 ) e 327I " 2 / A — wlogw . 



(2.10) 



In the regime of interest where A/(327r 2 ) is small, the dimensionless variables v,w and (p 



contain the very small factor e 



-32tt 2 /A 



times factors that are typically of order one. Hence, 



the physically interesting domain of values for v , w and ip is near the origin. Only for energies 
near the Landau ghost, v,w and ip become of order one. 

Hence, in the infinite N limit, the ('pion') field $ a (1 < a < N) is a free field with mass 
squared equal to M 2 . 

On the other hand, the composite field $ 2 (x) has as inverse propagator in euclidean 
space time given by 



A(ife) = A(0) 



(4vr) 2 Jo 



dx log 



(2.11) 
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where 



2 1 / \m 



2 



A<0) = -x-T^ los {7Jp) (2 ' 12) 

The relevance and properties of this propagator will be investigated below for the cases 
of unbroken m 2 > and broken m 2 < symmetry. 

It is clear from the gap equation ( |2.5| ) and fig. (|I|) , that 2 as a function of M 2 has a 
maximum at 

M 2 = M 2 = |m| 2 exp - l^j (2.13) 

for real M 2 . Notice that such mass scale is of the order of the Landau ghost. That means 

that it is outside the expected domain of validity of the ($ 2 ) , which can be used as an 
effective theory for energies well below the Landau ghost. 
At the maximum of 2 = 2 (M 2 ), we have 

*N*W> = ^ - 2j f ■ (2.14) 

This value is also of the order of the Landau ghost for our choice of m 2 . 

Therefore, for scales much lower than the Landau ghost, the (<3> 2 ) 2 theory is a sensible 
low energy effective description. In particular in our studies of the dynamics the scale of the 
expectaction value is typically of order 2 « \m\ 2 /X << 0q which is well within the regime 
of validity of the effective theory for weak coupling. 

As observed in pl,22|, the inverse function M 2 = M 2 (0 2 ) has two branches as a function 



of (J) 2 . However, the second branch, M 2 > Mq, (branch II in [|2T|) corresponds to pion masses 
beyond the Landau ghost. A consistent low energy effective theory emerges for energies well 
below Mq [See ref P2[, sec. IIIA]. This low energy condition requires that this scalar field 



theory be studied with a finite but large ultraviolet cutoff \m\ « A << \m\ e 167r / A and a 
positive coupling A. 



A. Analysis of the gap equation and the effective potential for m 2 > 

Equation fl2.5p and fig. 1 reveal that for m 2 > 0, the gap equation has a real and positive 
solution for < <p 2 < <p%, M 2 < M 2 < M 2 where M 2 corresponds to <p 2 = 0. The field <p 
becomes complex and hence unphysical for M 2 < M\. 

In fig. 2, the point 



2 



corresponds to the intersection of the lower curve (unbroken symmetry) with the (p — axis. 
W\ is thus the positive solution of the trascendental equation 

32tt 2 _32^ 

Wi log Wi = — e ^ . 

A 



[cfr. eq. (|2.10|) 1 which defines the value of M 2 . 



S 



Clearly, the effective potential is minimal for 2 = 0, M 2 = M 2 . This is the O(N) 
symmetric ground state with pions of mass Mf > 

The study of the field theory dynamics in this unbroken symmetry case for initial excited 
states [Tl[] shows that the (time dependent) effective mass tends to a nonzero value for late 



times. The pions are thus massive and the zero mode slowly dissipates its energy (there 
is a non-zero threshold for particle production) and attains the ground state (<f) = 0) very 
slowly, with an oscillatory behavior modulated in amplitude by a power law with anomalous 
relaxation exponent for asymptotically large time [[TT|j . 

We see that the static analysis through the effective potential agrees with the long time 
dynamics in the sense that the expectation value of the scalar field relaxes to the minimum 
of the effective potential, i.e. = with oscillations that reflect a non-zero mass. In 
addition, the effective pion mass obtained asymptotically in the dynamical situation depends 
on the initial energy density (which is conserved) [|TTJ. This dependence of course cannot be 



captured by a static analysis which only applies to the ground state. 

The <3> 2 propagator given by eqs. (|2.11( - |2.12|) has a tachyonic pole at A(m T ) = 0, if 
A(0) < 0. This occurs for M 2 < m 2 exp f^y- — 2^). That is, a tachyon is present unless 
the pion mass is of the order (or larger) than the Landau ghost. In particular, tachyons 




are absent if one chooses the second branch of the function M 2 = M 2 (0 2 ) [21]. As stated 
before, this scalar theory should only be interpreted to be valid as a low energy effective 
theory with log ^ = 0(1). The tachyon is then present but at a very large mass provided 
A is small. We find from eqs.( p.ll| - |2.12| ), 

A(k) k =°° - 

Thus the tachyon mass 

rriT = Tn exp 

\ A 

is beyond the range of validity of the effective (J& 2 ^) theory. 

In particular in our studies of the dynamics the typical expectation values are of order 
(f) 2 m m 2 /X and for these energy scales 0(m/\/~\) the ultra-heavy tachyonic states are 
never excited and they decouple from the low energy theory. We can thus safely ignore the 
existence of the tachyon. 




B. Analysis of the gap equation and the effective potential for m 2 < 

In this case, it is clear from eq. (|2.5| ) and fig. 2 that <fi 2 > 2 ^ for M 2 real in the interval 



A 

< M 2 < Mq with Mq given by eq.( |2.13| ). Equation ( |2.6p and fig. 2 show that the minimum 
of the effective potential corresponds to M 2 = 0, 2 = 2 -H-. That is, 



V 2 {w) > <^ 2 (0) 



A 

327T 2 



A 



for real < w < 1/e. The minimum of V e ff is at the point w = in the upper curve (broken 
symmetry) in fig. 2. 
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This yields a broken 0(7V)-symmetry ground state. As expected, we have N —1 massless 



pions which are the Goldstone bosons of the O(N) — > 0(N — 1) breaking [22] as can be seen 
from eq. ( |2.6| ) and fig. 2 which show that M 2 = is the minimum of the effective potential for 
M 2 > and below the Landau ghost. However, this value corresponds to a non-negative </> 2 
only for m 2 < 0. The fact that the minimum M 2 = of the effective potential happens in the 
real configuration space (for m 2 < 0) is in agreement with the asymptotic late time behavior 
found in refs. [jfl^MTIJ . There it was found that the dynamical (time dependent) effective 
mass of the pion fields tends asymptotically to zero for m 2 < (and initial energy below 
the false vacuum but above the ground state). The presence of massless particles allows the 
zero mode energy to dissipate efficently into them and the zero mode rapidly reaches the 
ground state. This happens very slowly in the m 2 > case (unbroken symmetry). In this 
case the pions are massive and energy dissipation occurs via non-linear resonances |TTJ . 

Now we come to one of the important features of the effective potential in the broken 
symmetry phase: it follows from eqs.( [2.9| ), ( |2.10| ) and ( |2.14j ) that there is no real solution 
M 2 (or w) for 2 > 0q. Therefore, M 2 and w as well as the effective potential become 
complex despite of the fact that the field (ft 2 is real. 

Whereas in a loop expansion this imaginary part is typically ignored with the hope 
that higher order corrections (loops) could remedy this situation, in the large N limit this 
argument does not hold. The physics associated with the imaginary part of the effective 
potential has been clarified in references [p~5| , [3~6]| . 

In the region in which the effective potential computed from eq. ([2.6|) or eq.( |2.9| ) turns 
to be complex, referred to as the spinodal region, the system is in a mixed state in which 
several phases coexist. The true equilibrium description of the system in the spinodal region 
is in terms of a Maxwell constructed effective potential, which is flat in the coexistence 



region < < — j- |T7|]. Although such construction is rather arbitrary at the level of 



phenomenological statistical mechanics (the lever rule), a consistent renormalization group 
analysis has revealed that the infrared limit of a carefully coarse grained effective action 
yields the Maxwell constructed effective potential |L7|. We also note that in this case of 
continuous symmetry the spinodal region extends all the way from the maximum to the 
minimum of the potential without a region of metastability. This must be contrasted with 
the case of discrete symmetry, wherein the spinodal ends in the limit of metastability before 
the minimum of the potential. The fact that the spinodal region reaches from the maximum 
to the minimum of the potential is a consequence of the presence of Goldstone bosons. 

At this point we summarize the relevant features of the effective potential in the large 
N limit to compare to the results obtained from the dynamical evolution 

• In the unbroken symmetry case with positive renormalized mass squared (m 2 > 0) the 
static effective potential has a unique minimum for a nonzero value of the squared mass 
M 2 of the excitations and vanishing expectation value of the field <fi. The dynamics 
agrees with a vanishing expectation value of the field reached asymptotically at late 
times with power laws [ ll| . In addition, the asymptotic effective mass obtained from 



the dynamical evolution depends on the initial expectation value of the field pT 



The analysis of ref. f22[ shows that the O(N) symmetry is broken in the effective 
theory for m 2 < where again, by effective theory we imply describing energies well 



below the Landau ghost scale Mq [see eq. fl2.13|) 1. 
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When the renormalized mass m 2 < 0, the minimum of the static effective potential is 
found at M 2 = and for a nonzero value of the field 0. The O(N) symmetry is then 
spontaneously broken and the massless scalars are the corresponding Goldstone bosons. 
The static effective potential computed from eq. (|2.6| ) or eq. (|2.9| ) is complex in the 
interval < <p 2 < 2|m| 2 /A; this is the spinodal region that signals phase coexistence. 
The spinodal region reaches from the maximum to the minimum of the potential unlike 
in the discrete symmetry case. 

The physics inside the spinodal region cannot be captured via a static description. An 
equilibrium description of this region can be given in terms of the Maxwell construction 
which leads to a flat effective potential for < <p 2 < 2m 2 / A. Such a description emerges 
in quantum field theory with a consistent treatment using the exact renormalization 
group [|T7|-|TT|. In the next sections we will study the dynamical Maxwell construction 
that leads to a flat dynamical effective potential in the spinodal region via the non- 
equilibrium dynamics. 

In section IV we study the dynamics of evolution in a broken symmetry case but when 
the initial value of is above the maximum of the tree level potential. In this case the 
dynamics is similar to that of the unbroken symmetry case. The effective potential is 
simply unable to describe correctly this physical situation. 

• Furthermore, the effective potential clearly illuminates the regime of validity of the 
scalar theory as an effective, cutoff low energy theory below the Landau ghost. Ob- 
viously, the physics is radically different when one considers pions and/or field expec- 
tations values at the energy scale of the Landau ghost. Our study of the dynamics 
in both the broken and unbroken symmetry phases will be restricted to energy scales 
~ \m\/y/\ and energy densities ~ m 4 /A i.e. well below the Landau ghost scale for 
weak coupling. Therefore this scalar theory is a reliable cutoff low energy model to 
study the non-perturbative dynamics for energy scales below the Landau ghost. 



III. NON-EQUILIBRIUM DYNAMICS IN THE LARGE N LIMIT: THE 

DENSITY MATRIX 



The field theory tool to study non-equilibrim dynamics is the quantum density matrix 
p{t) which obeys the functional Liouville equation 

^f = [#(*), P] 

where we admit a Hamiltonian that can depend explicitly on time as it would be the case 
in expanding cosmologies. A particularly illuminating representation for the field theoret- 
ical density matrix is the Schrodinger representation |14]|§, in which the field operator is 
diagonal. In such a representation 

pM-)^(-);t} = ^\p(t)m 

and the diagonal density matrix elements p[$(-), $(•); t] = V(<$>; t) give the functional prob- 
ability density for finding a field configuration with profile at time t. Therefore, this 
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representation allows to ask questions about the formation of semiclassical configurations in 
terms of probability functionals. This is the correct interpretation for the question of finding 
any arbitrary field configuration in the ensemble described by the non-equilibrium density 
matrix. 

We will study the general case of a biased initial condition leading to symmetry breaking 
dynamics to illuminate various features of the non-equilibrium dynamics and the implemen- 
tation of Goldstone's theorem. Therefore we write 



$(£, t) = (yN<f>(t) + a(x, t), ff 
($ a (f,t)) = y/N(f>(t) 5 11 ' 1 = Tr [p(t)$ a (x)] 



— * 

and \I/ is an 0(N — 1) vector. 

The canonical momenta conjugate to a and \& a are given by 

U a {x,t) = a(x,t) ; U%(x,t) = ^ a (x,t) 

and the Hamiltonian becomes 

H(t) = J d 3 x [\ul + + 1(V$) 2 + ^(W) 2 + V(<f>, a, *)} + ^ 2 (t) 

In the Schrodinger representation (at an arbitrary fixed time t Q ), the canonical momenta 
are represented as functional derivatives 

UJx) = -i 



UUx) 



5a(x) 
5 



and the Liouville equation for the density matrix becomes a functional differential equation 



m 



% < V,a\p(t)\a,V>= \ H[n a ,%,a,$;t} - H[U a ,U 9 ,a,^;tU < y,a\p(t)\a,V > . (3.1) 

The large N limit can be implemented by introducing an auxiliary field or alternatively 
using the following factorization 

(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 

where the constant term in ( |3.2|) is irrelevant as it will give a subleading correction in the 
large N limit. In this approximation, the potential ( [2.1Q becomes 



a 4 - 


-> 6(a 2 )a 2 + constant 


a 3 - 


-> 3(a 2 )a 


(* • ^) 2 - 


-> 2(^ 2 )^ 2 - (* 2 ) 2 


vj>V - 
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V{<f>, a, *) = V(VN<P) + a V x (t) + 



MM j2 , M 2 (t) iff2 



v x (t) 



N A / 2 2m r 



N(f)(t) 



A 



™ 2 + ^ 2 W + ^<* 2 >(*) 



A 



M\t) = m* + -<P\t) + — (¥')(*) 
M CT 2 (t) = m 2 + ^ 2 (t) + A ( ^ )(t) 



(3.7) 

(3.8) 

(3.9) 
(3.10) 



Since in this approximation the evolution Hamiltonian is quadratic (at the expense of 
a self-consistency condition), and the different O(N) components evolve independently, the 
density matrix factorizes as [[yj 



p(t) = p a (t) ® p*(t) 



(3.11) 



Furthermore, the density matrices can be chosen to be Gaussian, and will remain Gaus- 
sian under time evolution. 

Using spatial translational invariance, we can decompose the fields into their spatial 
Fourier modes and use these as a basis for the Schrodinger representation density matrix: 



a(x, t) 



d 3 k 



(Jr(t) exp(—ik ■ x) 



( 27r )3/2 "fc» 

$z(t) exp(—ik ■ x) . 



(3.12) 
(3.13) 



(2tt) 3 / 2 

The Gaussian density matrices acquire a simple form in terms of these Fourier modes 







p a [a,a,t] = Y[Af a ,k(t) exp j- df^ a k {t)u- k {t) o k (t)a_ k (t) 

k 

- B^ k (t)o k (t)(T- k (t) +ip ajk (t) [a- k (t) - cr-fc(t)]} , 
P*$, *] = IiM %k (t) exp L^lM $ k . $_ k _ # fe . 



B %k {t) $ k ■ 



(3.14) 



Only the component along the broken symmetry direction, o has an expectation value 
of the momentum, represented by p a , k (t) in ( |3.14j ), whereas the transverse directions have 
vanishing expectation value of the momentum. The kernels Bq, are scalars under 
0(N — 1). Although we will ultimately work at zero temperature, we have written the most 
general form of the density matrix, including a mixing term given by the kernels B a ^ jk 
above to display the most general form of the density matrix. When the mixing kernels 
-B<7,#,fc vanish, the density matrix collapses to a product of a wave functional and its complex 



13 



conjugate, i.e. it describes a pure state. Also note that hermiticity of the density matrix 
requires that the mixing kernels B a ^ be real. 

Now the Liouville equation ( |3.1| ) can be solved for the time dependent kernels by compar- 
ing the different powers of a, \P on both sides of the equation flOD . We obtain the following 



set of evolution equations for the kernels flU 

. A^fc = 1 
K,k 2 



^ = l(^ k -A* k ) (3.15) 



iA a , k = A% k - B 2 a>k - ul >k (t) (3.16) 

iB a)k = B a;k (A a;k -At !k ) (3.17) 

< k it) = k 2 + Mlit) (3.18) 

4> + M 2 (t)(p = . (3.19) 

■ A,l '- A ' (3.20) 



iA %k = A 2 -B 2 -u; 2 (t) (3.21) 



^-B*,fc = -B*,fc — (3.22) 

c^(t) = k 2 + M 2 {t). (3.23) 

Writing the kernels A k for cr, \& in terms of real and imaginary components A k (t) = 
■A-Rk(t) + iA Ik {t) and using the reality of the B k , we find the following invariants of evolution 



B k (t) B k {t ) 



A Rk {t) A Rk (t ) 

and 

Af k (t) 



(3.24) 



A Rk {t) + B k {t) 



constant (3.25) 

The condition (|3.25|) reflects unitary time evolution of the density matrix, i.e. the nor- 
malization is time independent. 

At this point we notice that the dynamics of a decouples from that of the rest of the 
fields i.e. it does not contribute to the evolution of the order parameter or to that of the 
fluctuations of the ^ fields (transverse directions). This is of course a consequence of the 
large N limit in which the dynamics is solely determined by the N — 1 transverse degrees 
of freedom. 

At this point we trace out the density matrix of o and study the evolution solely in terms 
of the order parameter and the density matrix for the transverse directions. 



The equations for the kernels can be simplified by defining [14 



A Rk (t) = A Bk (t) coth[Q k ] (3.26) 

B k (t) = --^t (3.27) 
smh[6 fc J 
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with Q k an arbitrary real parameter, which for an initial density matrix in thermal equilib- 
rium will be identified with the usual finite temperature factor 0*. = uJk/T. The complex 
quantity 



A k {t) = A R , k {t) + iA Ih {t) 
obeys a simple Ricatti differential equation 

iA\(t) = Al(t)-u 2 k (t) 
which can be solved by introducing a complex field 



A k (t) 



<pl(t)+u 2 k (t) <p m h {t) = 0. 



(3.28) 
(3.29) 

(3.30) 
(3.31) 



The initial conditions on the mode functions are determined by the initial values of the real 
and imaginary parts of the kernel A k (t = t ). Since an arbitrary time independent phase 
can be removed from the wave functionals, the imaginary part of the kernel A k (t) at the 
initial time can be set to zero without loss of generality. Thus we parametrize the initial 
conditions on the mode functions <fk{t) in the form 



<Pk(tc 



V2Wk 



V?fc(*)l 



to 




(3.32) 



With this choice the initial density matrix for zero mixing kernels corresponds to the product 
of vacuum wavefunctionals of harmonic oscillators of frequencies Wk and their complex 
conjugates, or if the mixing is of the thermal form it corresponds to the density matrix of a 
thermal ensemble of harmonic oscillators of these frequencies. 
We find in the large N limit 



(^ 2 (t)) = N 



d k . / \ 1 9 i ® & 
-. — -r \(PkW\ c °th — 



(3.33) 



The mode functions <fik{t) form a complete set of states allowing us to expand the Heisenberg 
operator for the transverse fluctuations as 



d 3 k 



a k (p k {t) e ik -* + a\ <p* k (t) e 



-ik-x 



(3.34) 



where the annihilation and creation operators will allow us to make contact with the Fock 
representation. With the initial conditions on the mode functions ( |3.32| ) we see that at 
the initial time to, to be chosen to — in what follows, the operators a k create quanta of 
frequency Wk out of the Fock vacuum state. 

It is convenient and illuminating to introduce the adiabatic creation and annihilation 
operators a/~(t) ; a\(t) and expand the Heisenberg field in the form 



ty(x) 



d 3 k 
(2tt) 3 



ctk(t) 



-i J dt'u h (t') 



Ak-x 



2u k {t) 



A Jdt'ui k (t') 



Ak-x 



2u k {t) 



(3.35) 
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The adiabatic operators diagonalize the instantaneous Hamiltonian and will provide a clear 
description of the particle production contribution to the total energy. The adiabatic op- 



erators are related by a Bogoliubov transformation to the operators a k 



,t 



that create 



excitations with respect to the Fock basis at the initial time. These operators are available 
whenever the time dependent frequencies are real. 

We now collect the final set of equations that determine the full dynamics of the density 
matrix in the large N limit 



<Pk{t) + 

+ 



2 , 



A 



2 2N X v ;/ 



d 3 k 
(27" 



3 ,V3 fc (t)| 2 coth 



= 



(3.36) 
(3.37) 
(3.38) 



These equations of motion are the same as those obtained in the Heisenberg represen- 
tation in the large N limit. The renormalization aspects have been studied in references 



10-141 to which we refer the reader for details. 



Although we have obtained the equations of motion in the general case of a mixed 
density matrix, we will consider in what follows the situation of an initial pure state by 
setting Q k = 0. 



IV. DYNAMICS FOR BROKEN SYMMETRY POTENTIAL 

In this section we study in detail the dynamics in a broken symmetry potential m 2 < 
in two important cases: i) V(t] ) < V(0) corresponding to [770 1 < v2, ii) V(r/ ) > V(0) 
corresponding to \r]o\ > y/2, where we define 

9r?i 4 N 
A 

V( V ) = \( V 2 - l) 2 . (4.1) 

We will analyze both of these situations numerically below. Our results can be summarized 
as follows. 

In the first case, V(r/o) < V(0), the energy is such that dynamics corresponds to a 
spontaneously broken situation in the sense that if r] 7^ the dynamical evolution results 
in a non-zero asymptotic value for 77(00). In these cases the effective mass of the transverse 
fluctuations vanishes asymptotically, i.e. these fluctuations are Goldstone bosons |6| TlH. As 



it will become clear from the dynamics, in this situation the Maxwell construction is realized 
through the non-perturbative production of Goldstone bosons. 

In the second case, V(t]q) > V(0), the energy is large enough that under the dynamical 
evolution the expectation value samples the minima ergodically. Its asymptotic value os- 
cillates around 77 = with diminishing amplitude and with a non-zero asymptotic effective 
mass for the fluctuations despite the fact that m 2 < 0. In this case particle production is a 
consequence of parametric amplification. 
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A. Equations of Motion: Numerical Analysis 



Here we set up the form of the equations of motion for the expectation value and the 
relevant modes in a form amenable to numerical analysis. 

It proves convenient to introduce the following dimensionless quantities in terms of the 
renormalized mass and couplings ( |2.2|J2.3| ): 

A 



m 



r 



Iml 2 — 



2N 



t = 0)> 



R 



\m\t , q 



A 



k 

Iml 



a, 



Iml 



1 
I 



1 1 

+ 



<?£(r) = ^2 (**(x,t)) R -(*\x,0)) R 

2 m 



A 32tt 2 

-A, 2 (t) 
2|m| 

( s(o) = o ; 







A 



m\ ip k (t) 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



Here {$ 2 (x, t))n stands for the renormalized composite operator [see eq. (|4.9|) for an explicit 
expression] . 

In the case of broken symmetry ( m 2 < 0) the field equations in the N = oo limit become 

IHI: 



V (r) = 



+ q 2 + M 2 {r) 



U q (r) 



Here, 



M 2 (t) 



(4.6) 
(4.7) 

(4.8) 



plays the role of a (time dependent) renormalized effective mass squared and £(t) is given 
in terms of the mode functions U q (r) by 
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gY,{T)=9 I q 2 dql\U q (r) \ 2 - — 



2q 3 



-Vo+V 2 (r)+gE(r) 



(4.9) 



The choice of boundary conditions is subtle for broken symmetry and a distinction must 
be made for the cases in which the initial expectation value of the scalar field is such that 
T]q > 1 or % < 1- in the first case, the initial frequencies are all real and vacuum type initial 
conditions can be chosen. In the second case, there is a band of unstable modes and the 
initial conditions must be chosen differently. The resulting dynamics in weak coupling is not 
very sensitive to the different choices of initial conditions in the second case 0. We shall 
use the following initial conditions for the mode functions HPf: for rjo > 1, we choose 



U q (0) 



u q (o) 



(4.10) 
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while for \rjo\ < 1 

U q (0) 

a, 



Vo 



for q 2 < ql = 1 



vl 



q 2 -l + r] 2 for q 2 > q 2 u ; < i] 2 < 1 . 



(4.11) 
(4.12) 



To these we add initial conditions for the zero mode given by 

7/(0) =7fo , 7)(0) =0 

Although the notion of particle is not unique in a time dependent non-equilibrium sit- 
uation, a suitable definition can be given with respect to some particular pointer state. 
Two definitions of the particle number prove particularly important for the discussion that 
follows. The first is the particle number defined with respect to the initial Fock vacuum 
state. This takes the following form in terms of the dimensionless mode functions we find 
the number of particles of (dimensionless) momentum q per unit volume to be given by || 



N q (r) 



^q\ U q{ T >\ + ^ 



1 

2 



(4.13) 



For the stable modes, i.e. for those for which the adiabatic frequencies are real at all 
times we can also consider the adiabatic particle number. It is given by flH|,pf 



Nf(r) = (a\{T)a q {r)) 



(4.14) 



w«(r) \U q {r)f 



\U q {r)\ 



1 

2 



where oj q {r) = yq 2 + A4 2 (t). Both definitions coincide at r = 0. 

We choose a initial state with zero temperature so the the initial occupation number 
N q (0) vanishes. It is easy to generalize our results to the case of initial states with a thermal 
distribution. 



V(?7o) < V(0), (|r?o| < V?*)' Evolution towards spontaneously broken symmetry 

states 



B. , 

states 

In this case the generic features of the evolution were studied in refs. RO]. If initially 
Vo — Vo — 0, then 77 remains at zero; this is a fixed point of the evolution equations. For 
early times, the mode functions grow exponentially as a consequence of spinodal instabilities 
and their amplitude become non-perturbatively large \U q (r)\ 2 w 1/g for r > t s w In \ j g (see 
eq.( |5.2| ) below and references P,pi]]). For rj 7^ we find that, whereas the effective mass of 
the fluctuations vanishes asymptotically leading to Goldstone bosons, the expectation value 
77 (r) tends asymptotically to a constant 77^ that depends on the initial conditions. Fig. 3 
displays this evolution for two different initial expectation values. It clearly shows that 77^ 
depends on the initial condition. Furthermore we can see from this figure that the effective 
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frequency of oscillation becomes smaller at larger times, consistent with a vanishing effective 
mass and the production of Goldstone bosons (see Fig. 5). 

The relation between the asymptotic value T]^ and the initial value r] Q is shown in fig. 4 
for < 7] < y/2. We find that vs. rjo can be fitted remarkably well (see fig. 4) by the 
power law behavior 



(4.15) 



with x a dynamical critical exponent. A numerical fit yields x = 0.25 . . . for small coupling. 
This exponent is rather insensitive to the coupling for weak coupling. Analogous dynamical 



exponents were found in ref. [II] for the late time evolution of the expectation value in the 
unbroken symmetry case. 



The law (|4.15|) is invariant under the exchange r/o <^> y 2 — rj 2 showing that ^oo(^o) is 
solely a function of the total initial (dimensionless) energy E = |(r/ 2 — l) 2 . In addition, we 
have ^(l) = 1 and /^(O) = 0, as expected, since these are fixed points of the dynamics for 
f](t = 0) = 0. Furthermore we find that the effective time dependent mass given by 



vanishes asymptotically as [JTT 



M 2 (t) 



A 



— cos 
r 



2r + 2 c In — 



7 



(4.16) 



with A ; c given in reference [llj|. Figure 5 shows Ai 2 (r). This asymptotic form of the 
effectiva mass translates into the following asymptotic behavior of r] (r) 



r)(r) 



1 + -M 2 (r) 
4 



+ 0(l/r 2 



(4.17) 



C. Dynamical Maxwell construction: 



The results of the numerical analysis lead to two important consequences that must be 
highlighted: 

• For |?7o| < \/2 all values of < < 1 can be reached asymptotically depending on 
the initial expectation value (energy). This is manifest in the expression ( 4.15 ) for the 
asymptotic expectation value and is the dynamical equivalent of the 'lever rule' for 
the Maxwell construction in equilibrium statistical mechanics. 

• Asymptotically the effective mass squared, i.e. the "restoring" force in the evolution 
of tj(t) vanishes as given by eq. (|4. 16|) (see fig. 5), compatible with massless Goldstone 
bosons. 

• The states for which < [7700] < 1 are excited states with large energy densities, 
the smaller 7700 the larger the energy density. Most of the energy is stored in long- 
wavelength condensates as a result of non-perturbative particle production via spinodal 
instabilities. 
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These features of the dynamics combined lead to the conclusion that all of the broken 
symmetry states with asymptotic equilibrium values < < 1 are allowed and that 
the dynamical potential associated with these states is flat since the "restoring force" for 
the evolution of the order parameter is Ai 2 (r) which vanishes asymptotically as a result of 
Goldstone's theorem. Therefore we see the emergence of a dynamical Maxwell construction. 
The energy initially stored in one mode, (the expectation value or zero mode) is transferred 
to all of the other modes but mainly to those that are spinodally unstable. The amplitude 
of these modes becomes non-perturbatively large ~ 1/g. 

This should be contrasted with the results from the static effective potential that predicts 
that the only equilibrium broken symmetry solution which is physically acceptable is \r]\ = 1, 
states with < 1 770 1 < 1 lead to an imaginary part, which as described previously is the 
hallmark of the spinodal instabilities leading to non-perturbative particle production. 

In fact, we can make contact with the static effective potential picture in the following 
way. 

Since we are considering translationally and rotationally invariant states, the expectation 
value of T^ v takes the perfect fluid form. The energy density, i.e. expectation value of the 
Hamiltonian (divided by the volume and N) in the broken symmetry case is given by || 



NV 
In ref. 



2 Ao 4 _ Ao. ( ^ 2) 2 
2 y 8 y 8N 



k 2 alk 



\ti(t)\ + "t(t)\Mt)\ 2 (4.1 



2 2 T ' 8 T 8N 1_ ' ' 4vr 2 
it was proven that after a constant subtraction oc A 4 the energy is finite , 



furthermore ( 4.18 ) is conservedhy use of the eqs.(4i;) and 

Using eq. Q4.14j ) and eq.( [4.9| ) we write the last term in ( [4.18| ) in the following form 



1 

in- 



k 2 dk 



\ip k (t)\ 2 + u 2 k (t)\Mt)\'' 



£u + 



e v = 



1 ^ k 2 dk 



2 Jo 2vr 2 



\y k {t)\ 2 + uj 2 k {t) |p fc (t)| 5 



(4.19) 
(4.20) 



where k 2 is the maximum unstable frequency in the broken symmetry case for the case 
r) 2 < 1. The term Ejj is the contribution to the energy from the unstable modes with negative 
squared frequencies, and N% d {t) is the adiabatic particle number given by eq. ( |4.14|) . We have 



written explicitly an upper momentum cutoff, so as to absorb the divergences in the mass 
and coupling renormalizations. Using the large iV renormalization condition |6|,|l~4|j 



2 , Ao , 2 . A ,^2\ 

rrir, H d> H (W } 

2 2N X 1 



m 2 + 



A 



+ —(&) R 

T 2N \* I R 



(4.21) 



and the renormalization prescription used for the effective potential given by eqns. ( |2.2| . 
|2.3|) we find the once subtracted (by a time and field independent term oc A 4 ) renormalized 
energy density divided by N to be given in the A — > 00 limit by 



+ 



^ 2 0\,2,n M 2 (t)\m\ 2 MHt) 

- + -7-M 2 (t) + £ F (t) 



A 



k u M 2 {t) 
32tt 2 



(k 2 u + M 2 (t) 



2A 

M 4 (t), Ik,-, 



M 2 {t) = -\m\ 2 + ^ 2 (t) + ^(&(x,t)) R 



16 



\m\ 2 \fe 



U V2 {k 2 u + M\t)- 
2 

(4.22) 
(4.23) 



20 



where M 2 (t) is the effective time dependent mass of the transverse fluctuations and 
1 f A 



10 



e F (t) = — \i k 2 dk \<p k (t)\ 2 + u 2 k (t) \<p k (t)\ 2 + / k 2 dkcu k (t)N^(t)\ (4.24) 



contains all of the dynamical information on the spinodal instabilities as well as of adiabatic 
particle production in the broken symmetry case. This is because the mode functions for 
the unstable modes will grow exponentially |16| , p7| during the early stages of the dynamics. 



The total energy is finite, conserved and real. 

We can now make contact with the static effective potential formalism. Recall that the 
static effective potential is a quantity defined for time independent expectation value of the 
field. For this purpose, we identify the saddle point solution of the Lagrange multiplier M 2 
with the effective mass of the transverse fluctuations, which in the dynamical case is M 2 (t) 
given by eq. ( [4.231 ). Comparing the effective potential eq. ( |2.4j ) with the above expression 



for the total energy, we see that we obtain the effective potential from the energy by: i) 
setting <f) = , ii) neglecting the contribution ep, iii) setting k u = 0. Clearly setting k u = 
leads to an imaginary part of the effective potential whenever the effective mass is negative. 
This is the region of spinodal instabilities in the case with the initial expectation value is 
such that ?7q < 2, and whose dynamics is described precisely by ep. Even in the case in 
which the initial condition is such that rj 2 , > 2, for which k u = 0, the effective potential 
completely misses the term ep which is now completely determined by the adiabatic particle 
production. This is the regime of the dynamics in which particle production is a result of 
parametric amplification rather than spinodal instabilities p|j6|Jn| . 

We thus see that particle production via spinodal instabilities for t)q < 2, which is 
accounted for in the contribution ep, solves the problem of the imaginary part and is re- 
sponsible for a dynamical Maxwell construction. For t)q > 2 in which case the initial state 
has an energy larger than the top of the potential barrier and the dynamics can probe er- 
godically the minima, parametric amplification leads to profuse particle production which 
is accounted for by the contribution ep now with k u = 0. Therefore in both cases when 
particle production is important, either in the spinodally unstable region or for large energy 
densities when particle production is via parametric amplification, the effective potential 
provides a misleading picture. However these are the important regions for the study of 
phase transitions or large energy densities. 



D. V(r/o) > V(0), (|ryo| > V2): symmetric evolution 

We now focus on the case of the broken symmetry potential, but when the initial expec- 
tation value is higher in the potential hill than the maximum, i.e. |?7o| > V2. In this case 
the total energy density is non-perturbatively larger than V(0). Classically the expectation 
value of the scalar field will undergo large amplitude oscillations between the two classical 
turning points thus probing both broken symmetry states ergodically. The classical evolu- 
tion is completely symmetric in the sense that averaging the evolution of the expectation 
value over time scales much larger than the oscillation period will yield a vanishing value, in 
contrast with the classical evolution for |??o| < \/2- This is not the situation that is envisaged 
in usual symmetry breaking scenarios. For broken symmetry situations there are no finite 
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energy field configurations that can sample both vacua. In the case under consideration with 
the zero mode of the scalar field with large amplitude and with an energy density larger than 
the top of the potential hill, there is enough energy in the system to sample both vacua. 

Parametric amplification transfers energy from the expectation value to the quantum 
fluctuations [p],|6|,|TT|| . Let us assume that a large fraction of the energy in the expectation 
value has been transferred to the quantum fluctuations such that the amplitude of the 
expectation value now is \t](t)\ < y/2. Since the quantum fluctuations react back on the 
dynamics of the expectation value and now there is an energy density oc 1/g stored in the 
quantum fluctuations, the quantum fluctuations act as a bath with large energy density 
that forces the expectation value to continue to sample ergodically both symmetry breaking 
vacuum manifolds. 

Fig. 6 displays minus the effective mass term, — A4 2 (r), as a function of time. Clearly 
the effective mass is positive, translating into the fact that X(^f 2 )/N > \m\ 2 as a consequence 
of particle production via parametric amplification. Fig. 7 displays 77 (r) for 7/(0) = 5 ; g — 
10~ 5 . 

This situation is qualitatively similar to that of equilibrium finite temperature field theory 
at temperatures larger than the critical in the sense that the mean root square fluctuation 
of the field is larger than the value of the field at the minima. Although qualitatively similar 
to finite temperature in equilibrium, the situation is quantitatively different because the 
distribution function for the quantum fluctuations is not thermal pjll|. A more proper 
description would be that the total system composed of the expectation value and the 
quantum fluctuations is described by a microcanonical ensemble at fixed energy. However 
treating the quantum fluctuations as a bath in interaction with the (one) degree of freedom 
described by the expectation value, one can think of an alternative canonical description 
in which the expectation value exchanges energy with the bath. With \(^> 2 )r/N > \m\ 2 
the fluctuations of the bath are large enough to allow the expectation value to sample both 
broken symmetry states with equal probability. This alternative canonical interpretation 
does not imply that the distribution functions are thermal since there is no restriction of 
maximization of any form of entropy. 

The evolution of the expectation value is damped because of the transfer of energy to 
all the other modes via the mean field. In the present dynamical case, this 'symmetry 
restoration' is just a consequence of the fact that there is a large energy density in the initial 
state, much larger than the top of the tree level potential, thus under the dynamical evolution 
the system samples both vacua equally. We note however that symmetric evolution would 
ensue even in the classical theory without quantum backreaction simply because the energy 
density is large enough. The initial state is not a broken symmetry state. It has a large 
enough energy density to allow the expectation value to overcome the potential barrier. 

Thus the criterion for symmetric evolution when the tree level potential allows for broken 
symmetry states is simply that the energy density in the initial state be larger than the top 
of the tree level potential. That is when the initial amplitude of the expectation value is 
such that V(i]o) > V(0). In this case the dynamics is similar to the unbroken symmetry case, 
the amplitude of the zero mode damps out, transferring energy to the quantum fluctuations 
via parametric amplification and non-linear resonances fTT||. We see that the amplitude of 



the zero mode tends asymptotically to zero very slowly in a manner similar to the unbroken 
symmetry case (m 2 > 0) p|JTT|] while it oscillates symmetrically around 77 = 0. 
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These points are clearly illustrated in figs. 6 and 7, M. 2 {t) and r/(r) for rjo = 5 > y/2 
[and hence V(t]q) > V(0)], where V[rj) is given by eq. ( |4.1| ) and g = 10~ 5 . We find the 
typical behaviour of unbroken symmetry P,|IT|. It is clear from the figure the presence of 
two frequencies: one corresponds to twice the asymptotic meson mass JA(oo). 

We see from fig. 6 that M. 2 {t) tends asymptotically to 



This value follows from the unbroken symmetry result [[11]] upon changing the sign of the 
tree level mass squared +1 — ► — 1. 

We emphasize again, that the fact that we find symmetric evolution for a negative tree 
level mass squared is not a surprising result, it is due to the high energy of the initial 
states considered ~ 0(1/ g) >> 1. Since the dynamical evolution sampled both vacua 
symmetrically from the beginning (even in the classical theory), there never was symmetry 
breaking in the first place, and 'symmetry restoration' is just the statement that the initial 
state has enough energy density so that the dynamics probes both vacua symmetrically 
despite the fact that the tree level potential allows for broken symmetry ground states. 

Energy conservation and mode mixing through the mean field leads to a damped evo- 
lution for the expectation value but non-perturbative particle production via parametric 
amplification, results in a highly excited quantum state in which the energy has been dis- 



tributed mainly throughout the states in the parametrically unstable band(s) PJTT|j. The 
sign of the tree level mass squared which is so crucial near the ground state becomes here 
almost irrelevant. 



V. QUANTUM PHASE ORDERING DYNAMICS 

Phase ordering refers to the process phase separation after a sudden cooling from the 
disordered, high temperature phase into the low temperature phase. During this process the 
system breaks up in ordered domains, this is the essence of Kibble's mechanism [p3|-|25l for 



the formation of defects after a phase transition. Whereas the dynamics of phase ordering 
is fairly well understood in condensed matter physics |30| , |3l| , a microscopic quantum field 



theoretical treatment is still missing. Analogies with condensed matter physics had been 
used to provide some dynamical estimates in quantum field theory p6fl , but a full non- 



equilibrium, dynamical treatment beginning from a microscopic quantum field theory is just 
beginning to emerge [p^-p9|. Analogies with the condensed matter situation led to the 



study of classical scaling solutions for the dynamics of texture collapse and relaxation in the 
non-linear sigma model in a matter-radiation dominated cosmology [^.Q • 

In this section we provide a description of the process of phase ordering, implementing 
the non-equilibrium formulation in the full quantum field theory (QFT). 

Consider the situation in which the system has been quenched from a symmetric (disor- 
dered), high temperature phase into the low temperature (ordered) broken symmetry phase. 

In QFT this quench can be introduced by considering an initial quantum mechanical wave 
function(al) prepared at t — centered on top of the potential, i.e. a wave function(al) that 
is localized near the origin in field space |15|j27|JT^] . This wave functional is then evolved 
in time with the Hamiltonian with a potential of the broken symmetry form. This wave 
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function(al) will evolve in time by spreading out to sample the equilibrium states, which 
results in a growth of correlations and fluctuations. It is precisely in this situation that the 
density matrix formulation and its probabilistic interpretation offer the most enlightening 
description. The early time dynamics relaxing the quench assumption has been been recently 
studied in reference p4| . 

In the large N limit the concept of topological defect is rather ambiguous, and as we will 
see, the important concept is the formation of long-wavelength, large amplitude (semiclas- 
sical) field configurations during the evolution. 

Let us consider the situation of a critical quench, with initial conditions of vanishing 
expectation value and first derivative at the initial time, i.e. r/(r = 0) = ; i)(t = 0) = 
0. Since this is a fixed point of the dynamics for the expectation value, it will remain 
zero all throughout the evolution. From the probability distribution given by eq.( |5.15| ) 
and the initial conditions on the mode functions ( }4.11|) -( }4.12|) , we see that at the initial 
time the most probable configurations are those whose spatial Fourier transform are of 
the order 1/^/1^. Therefore the probability for finding semiclassical long wavelength field 
configurations with non-perturbatively large amplitudes, i.e. oc 1/^fg, in the quantum 
ensemble is exponentially suppressed by the inverse of the coupling constant, i.e. P^sc] oc 

However, from the evolution equation for the mode functions ( [4.7| ), we see that the 
modes with wavevectors in the spinodally unstable band grow exponentially and their am- 
plitude becomes of order at the spinodal time r s ~ ln(l/g). Figures (11)-(12) show 
y/g\U q (T = 100,200)| vs. q respectively for g = 10~ 7 ; t](t = 0) = t](t = 0) = 0. We 
infer that for tau > t s there is a unsuppressed probability for finding long wavelength 
k < \m\ (q < 1) large amplitude oc 1/y/g field configurations in the ensemble. The interpreta- 
tion of this observation is that spinodal instabilities lead to the formation of non-perturbative 
semiclassical long-wavelength field configurations with amplitude ~ 1/ yfg in some functional 
direction. This is consistent with the fluctuations sampling the minima of the potential, 
that is (ty 2 )/N ~ \m\ 2 /X or mean square root fluctuation of the field ~ |m|/\/X 

An important concept in the dynamics of phase ordering is that of the time dependent 
correlation length. In condensed matter systems the process of phase ordering at late 
times (after initial transients shortly after the quench) is described in terms of correlation 
functions that show a scaling form in terms of a dynamical correlation length which at long 
times is independent of the (short) initial correlation length in the disordered phase. A 
scaling solution has also been found |32| in the classical large N limit but in matter or 
radiation Friedman- Robert son- Walker cosmology. 

We now turn our attention to a detailed study of the equal time correlation function to 
understand the emergence of scaling. 



A. Correlation Functions 



Consider the case of a critical quench i.e. ($) = ($) = , i.e. 77(0) = fj(0) = 0. In 
terms of the dimensionless coordinate distance r = \m\x and the dimensionless quantities 
introduced before (|4.2|44.5|) the equal time correlation function becomes 



< $ a (f,t)$ 6 (0,t) > = C(r,r) 8 



a.b 
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= H 2 S a ' b J ^ <** |^(r)| 2 = ^ jfffdg sin(gr) |^(r)| 2 . (5.1) 

with r = |f|. 

We display in fig. 10, grC(r,r) as a function of r for various values of r. This figure 
reveals several remarkable features that we analyze in detail below. 



B. Intermediate time: r < r. 



For times r earlier than the nonlinear or spinodal time t s we can neglect the back-reaction 
of the quantum fluctuations in the evolution equations for the mode functions. This is the 
linear regime and we can use the spinodal behavior of the mode functions in order to describe 
the behavior of C(r, r). When r « r s , the quantum fluctuations become non-perturbatively 
large and the back-reaction, encoded in the term <?S( r ) * n the evolution equations becomes 
comparable to the tree level term. For weak coupling, the spinodal time scale is given by 
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(5.2) 



and the mode functions for r < t s can be reliably approximated as |Tl 

1 



V / l-<7 2 -*(l + <7 2 ) 



r^/T- 



+ e 



for wavevectors in the spinodally unstable band < q 2 < 1. Inserting this expression into 



eq.(|5.1|) yields upon saddle point integration p7|| 

1 



C(r,r) 



1 + C 



r 



(5.3) 



16 (tt r) 3 / 2 

Eq.( |5.3|) is in remarkable agreement with our numerical calculations in its range of validity 
|TT| . The exponential decay with r in eq. (|5.3j ) can be interpreted as the emergence of a 
dynamical correlation length 

£ e = 2v^ (5.4) 
The correlation function in this regime can be written in the form 



C(r,r)=g(r)F 



T 



(5.5) 



which, however, is not a true scaling form, since the function g{r) is not a power of the 
dynamical correlation length. This behavior is reminiscent of a diffusive process and similar 
to the results of phase ordering kinetics for non-conserved order parameter in the non- 
relativistic case [p0|1 , which will be briefly reviewed in the next section. The 'diffusive' type 
behavior in this case is understood from the fact that the modes that grow the fastest are 
those for which q « 1 for which a non-relativistic approximation of the modes is valid. 

As we will see below, however, the \fr behavior of the dynamical correlation length is 
only valid for r < r s ; there is a crossover to a linear scaling law and a true scaling solution 
at longer times. 
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C. Short Distance Analysis 



The correlator for small distances can be obtained by expanding sin(gr) in eq.( |5.1|) and 
restricting the integration to the band of spinodally unstable momenta q < 1, since the 
contributions from q > 1 are perturbatively small corrections. Doing this yields: 



C(r, r 



47T S 



„2 rl 



^(r)--^ qUq \U q (r)\* + O [r 



For late times we can use the sum rules arising from energy conservation |JTT 

0E(oo) = l , f l qUq \U q (oo)\ 2 = -jL . 

Jo Ag 

Therefore, 



C(r,r 



r^0,T>>l 



47T 2 5( 
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+ 0(r 



This behaviour is in perfect agreement with our numerical results for short distances. 



D. Large distance and late time behavior 

For r >> 1 the integral is dominated by the small q-region, and since the only important 
modes are those in the spinodally unstable band q < 1, changing variables to Q = qr the 
correlation function can be expressed as 



rC(r, 



T 



2vr 2 r 2 



QdQ sinQ |C/ fi (r)p 



(5.6) 



For r >> 1, only the mode functions U q (r) with very small q are relevant. Their long 
time behavior can be obtained from the following reasoning. At late times r >> r s the 
effective mass term vanishes asymptotically as given by eq. ( 4.16|) Jll]. 

Since the effective mass vanishes asymptotically at long time as given by ( |4.16| ) [ H | there 
are no secular terms in the perturbative solution for the q = mode |ll] so that the mode 
with q = must behave as (Tl| 



Uq{t) = L + Kt + small oscillatory terms (5-7) 

The Wronskian guarantees that neither of the complex coefficients, L, K can vanish ||11|| . 
Figs.(ll)-(12) show both the real (y/g Ur j0 (t)) and imaginary (yfg Ui tG {r)) parts of y/g U (r) 
respectively for 77(0) = fj(0) = 0; g = 10~ 7 where it can be seen that L ~ K ~ Oi\j ^fg). 
This linear growth with time signals the onset of a novel form of Bose Einstein condensation, 
in the sense that at very large times the zero momentum mode will become macroscopically 
occupied |lTJ. Since the mass term vanishes asymptotically, the solutions to the equations 
of motion will be of the form 



U q (r) 



A q e 



ll[T 



+ B e 



-iqr 



The equation of motion for the modes is analytic in the momentum variable q, therefore 
analyticity in q fixes the small q behavior at long times to be of the form 
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K 

Uq«i\T) = L cos(gr) + — sin(gr). (5.8) 

Figs. (8)-(9) show the behavior of the mode functions at two large times, for very small q, 
the behavior of the modes as a function of q is very well described by the form ( |5.8| ) with a 
numerical error in the small q region less than a few percent. 

At large time and distance (r; r >> 1) the correlation function ( |5.6| ) becomes 

n( v r»l,T-oo |^| 2 f°°dQ . . 2 ( n T\ \K\ 2 

rC(r,r) = -ptt -77 sin Q sin [Q - = — — 6(2r - r) . (5.9) 
2tt z Jo Q \ r / qtt 

The correlation function falls off as 1/r for r < 2r and vanishes for r > 2r. which is 
precisely the behavior found numerically and shown in detail in fig. 10 which reveals that 
the correlation function vanishes for 

r>2(r-T s ) + c (5.10) 

with t s the spinodal time and c ~ 2 for small coupling. A similar analysis using the mode 
functions (|5.8|) leads to the cos(2r)/r behavior of A4 2 (t) at large times. 

The 1/r fall-off can be interpreted as a correlation function of massless Goldstone modes. 
However, these are not free modes, as their correlator would have a 1/r 2 fall off. 

The interpretation is the following, for times r > t s there is a zero momentum condensate 
formed by Goldstone bosons travelling at the speed of light and back-to-back. That is, 
massless particles emitted from the points (0, r) and (r, r) form propagating fronts which at 
time r are at a distance r — t s from the origin and from r, respectively. 

These space-time points are causally connected for 2(r — t s ) > r. Otherwise, the corre- 
lator vanishes. Causality allows the addition of a positive constant c as in eq.( 5.10| ). 



The root mean square fluctuation of the field inside this propagating front is non- 
perturbatively large ~ \m\/ ' ^fg reflecting the process of phase ordering. The fact that 
the root mean square fluctuation of the field inside these domains means that although 
the expectation value of the order parameter vanishes, inside these growing domains the 
fluctuations sample the broken symmetry minima. These then are ordered domains with an 
average expectation value inside each domain to be given the mean square root fluctuation 
and therefore saturated at the equilibrium values at the minima of the potential. 

For r < 2(r — r s ) + c we find numerically that 

C(r,r)^— (5.11) 
9 r 

for 77(0) = 77(0) = and small coupling. This result is in complete agreement with the 
estimates ( |5.9| ) since from figures (11)-(12) (linear growth) we infer that g\K\ 2 pa 2. 

Inserting the low q behaviour of the mode functions (|5.8|) into eq. (|4.9|) for yields for late 
times 

S(r) r =°°S(oo) + /H!l27 "' 



Comparing this with the results of ref. ]TT], we see that this approximation yields the cor- 
rect asymptotic behaviour for S(r) except for the logarithmic phases in the 1/r term [see 
eq. (P^|) 1. 
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The asymptotic form of the equal time two-point correlation function given by eqn. ( |5.9| ) 
reveals a true scaling solution. Introducing the dynamical correlation length £(r) ~ r and 
defining the variable x = r / £(t) it can be written in the form 

C(r,r) = [ar)]- 2{1 - z) F[x} (5.12) 

with the anomalous dynamical exponent z = 1/2 (the naive scaling length dimension of the 
field is 1) and the scaling function 

Fix] oc -9(2 - x) (5.13) 

x 

This scaling solution reveals the non-perturbative buildup of a large anomalous dimen- 
sion. 



E. Correlation length 

We can compute the correlation length at time r using the standard expression, 

= Jd 3 r \r\ C(r,r) 
^ [T) ~ fd 3 rC(r,r) 

There are two different time regimes: i) the early-intermediate time r < t s and ii) the 
late time r >> r s with different behavior for the correlation length. Thus there is a crossover 
in the scaling behavior of the correlation function at the spinodal time scale t s , with the 
correlation length behaving as 

£(r)oc{ ^ ;° rT ^ (5.14) 
I t — T s i or r > t s 

For times larger than the spinodal time, when the spinodal fluctuations have been shut- 
off by the growth of the self-consistent field, and the dynamics is driven by the non-linear 
evolution, the correlation length grows monotonically with time becoming infinite as t — > oo; 
this is a consequence of the presence of massless particles. Moreover, the vanishing of the 
correlation function for r > 2r is consistent with causality and tells us that the correlation 
length must be of the order of magnitude given by eq.( |5.14|) for at long times. 

Only for infinite time the correlation length takes the value infinite as one would expect 
in the presence of massless particles. The time effectively acts as a infrared cutoff making £ 
finite for finite time. 



F. Where are the Defects? 

As mentioned above, there are no defects as such in the O(N) model in the large N limit. 
However, what we are seeing is that this is not necessarily the relevant question, in terms of 
the phase ordering. What we should look for are long wavelength, non-perturbatively large 



configurations [12 
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These are always present in the ensemble described by the density matrix of the sys- 
tem. However, if our system were just a thermal one, such configurations would be highly 
suppressed, since their energies would be proportional to 1/A. 

What we see in our evolution, though, is that spinodal instabilities drive the growth of 
long wavelength, non-perturbatively large configurations. We should then expect to find 
them with unsuppressed probability in the ensemble after the spinodal time, say. A way to 
quantify this is to use the probability distribution obtained from the density matrix. 

The probability of finding a particular field configuration with spatial Fourier transform 

at time t in the ensemble is given by 

'^i^^rsij (5 - i5) 

which shows clearly that at a given time t, the most likely configurations to be found in the 
ensemble are those whose spatial Fourier transform are given by (pk(t). 

The mode functions </?&(£) become non-perturbatively large, oc 1/y/X because of spin- 
odal unstabilities. Therefore even when at the initial time field configurations with non- 
perturbatively large amplitudes are exponentially suppressed for weak couplings, at times 
of t s oc ln(l/A) large amplitude field configurations will be unsuppressed in the quantum 
ensemble. 

This interpretation is further supported by the results of reference [EL 



G. A classical but stochastic description 

At the end of the stage of linear instabilities at r « t s we have seen that the mode 
functions in the spinodally unstable band < q < 1 attain non-perturbative large ampli- 
tudes ~ 1/ y/g, whereas the larger wavelength modes remain perturbatively small. All of the 
correlation functions are dominated by the modes that achieved non-perturbatively large 
amplitudes through spinodal decomposition. This suggests that the field should be decom- 
posed as the sum of a classical part corresponding to the amplified modes and a quantum 
mechanical part corresponding to those modes with perturbative amplitudes 

= $ ae (x,t)+$(x,t) 

*-<*• () = jL<m wh ft ^ Jt) + Ai (5 - 16) 

d 3 k 



(2tt) 3 / 2 



where fj: sc are the 0(1/ y/g) parts of the mode functions, whereas (pj: qrn are the 0(1), 
i.e. perturbatively small part. Clearly the split between the semiclassical and quantum 
components of the field is ambiguous. However, this ambiguity is of O(g) corresponding to 
extracting the small components of the mode functions and therefore remains perturbatively 
small. 
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Whereas the a are quantum mechanical operators with usual commutation relations, 
the variables A are classical. However in order to reproduce the non-perturbative contribu- 
tions of the correlation functions these variables must be described in terms of a Gaussian 
stochastic probability distribution so that 

<< >>=<< A~ »= 

«^4»=%, (5.17) 

The stochasticity of these semiclassical variables is necessary to reproduce all of the correla- 
tion functions that are obtained from averages with the quantum mechanical density matrix 
given by ( |3.14| ). Because in the large N limit the quantum density matrix is Gaussian, the 
semiclassical probability distribution function is also Gaussian. 

We see that a stochastic description arises naturally from the time-evolution of the 
system, and is not put in by hand. 



H. Summary of quantum phase ordering dynamics 

We summarize here the results from this section to provide a picture of the process of 
quantum phase ordering in the large iV limit for the case of a critical quench, i.e. 77(0) = 
t)(0) = 0. 

• The process begins with the exponential growth of spinodally unstable modes. The 
equal space-time two point correlation function grows and begins to compete with the 
tree level terms in the equations of motion. This is the region of linear instabilities 
that lasts up to the spinodal time t s ~ l/\m\ ln(l/g) when g"S(r) ~ 1. The correlation 
function is of the scaling form, with a time dependent correlation length that scales 
as oc \/t. The growth of quantum fluctuations are interpreted as the formation 
of semiclassical coherent condensates with wavelengths k < \m\, with unsuppressed 
probability of being represented in the quantum ensemble. 

• At t ~ t s the quantum backreaction of the modes produces a crossover and for t > 



t s the evolution is non-linear ||11|| . The effective mass of the fluctuations vanishes 
and the distribution functions reveal non-perturbative long-wavelength condensates 
of Goldstone bosons. A novel form of Bose Einstein condensation is manifest in the 
linear growth with time of the zero momentum mode function. This soft condensate 
of Goldstone bosons in turn leads to a correlation function that falls off as 1/r and is 
cut-off by causality at r > 2(t — t s ). We interpret this result as that the process of 
phase ordering occurs via the formation of a domain or bubble that grows at the speed 
of light for t > t s inside which there is a non-perturbative, semiclassical condensate 
of Goldstone bosons. Inside these domains the mean square root fluctuations are 
~ \m\/\f\ and therefore sampling the equilibrium minima of the potential, i.e. the 
order parameter averaged over this domain saturates at the equilibrium value. The 
equal time correlation function falls off as 1/r in contrast to free field behavior 1/r 2 . 

For t >t s the non-perturbative aspects of the dynamics can be described by a semi- 
classical but stochastic distribution function. The semiclassical description is achieved 
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through the field expansion in terms of the modes in the spinodally unstable band with 
coefficients that are interpreted as stochastic variables with a Gaussian probability dis- 
tribution function. This description reproduces the non-perturbative contributions to 
all of the correlation functions in the large N limit. 



VI. COMPARISON WITH CLASSICAL PHASE ORDERING KINETICS: 

We now compare the results obtained above on quantum phase ordering dynamics to 
those previously found in the literature but in different contexts but within classical stochas- 
tic field theory: i) Non-relativistic condensed matter: phase ordering kinetics for non- 
conserved order parameter in the large N limit [30|H]|, ii) Cosmology: the relaxation of 



texture-like configurations in a Friedmann- Robert son- Walker background metric in the non- 



linear sigma model in the large N limit |j32|| . Our main purpose is to compare and contrast 
our results to those obtained within the classical large N approximation and to point out 
that the body of results found here describe a fundamentally different set of phenomena and 
scaling. 



A. Phase ordering kinetics in condensed matter: 

The process of ordering after a quench in condensed matter physics has been studied 
both theoretically and experimentally very intensely during the last decade (for a thorough 
review see ||30|| ). Amongst the theoretical approaches, the large N limit stands out as one of 
the few that yields to an exactly solvable description of the dynamics and used as a yardstick 
against which to compare other approaches [E(J. Let us now briefly summarize the most 
relevant features of the solution in the large N limit. For non-conserved order parameter, 
the process of phase ordering is studied by means of the time dependent Ginsburg-Landau 
equation (TDGL) 

9^,t) = r 5F 
dt 6${x,t) 

Here T a phenomeno logical dissipative coefficient, which in this model can be absorbed in a 
definition of time, and the Ginzburg-Landau free energy F is given by 



d x 



Uv$(x,t)f + V{${x,t)) 



(6.2) 



with V(4>) given by eq. ( p. ID (with m 2 , < for broken symmetry). The TDGL equation is 
purely dissipative and a simple calculation of the energy, which is given by total spatial 
integral of the free energy density, shows that it is a monotonically decreasing function of 
time. 

In terms of the rescaled field <fi 2 = (2iV|mo| 2 /Ao) ff and rescaled time and spatial variables 
|mo| 2 rt = r ; |too|£ = z the TDGL equation now reads 

^ = V 2 v + v~(v) 2 V (6.3) 
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The large N limit is implemented by a Hartree-like factorization of the non-linear term 
HI V 2 ~ > (tf 2 ) where the average is over the initial probability distribution function. The 



initial state is assumed to be completely disordered and described by a Gaussian distribution 
function with 

(rf(z,0)) = ; (f|"(z,0)i)^,0)) = A^(5 3 (?-?) (6.4) 



with A a constant determined by the high temperature short ranged correlations [pOjpl 
This factorization leads to a simple equation 

M\t) = 1 - (ff) (6.5) 

but with a self-consistency condition much in the same manner as the large N description 
in quantum field theory, notice that the effective mass term _M 2 (r) is similar to that found 
in eq. (O). However, notice that eqs.( |6.5| ) are not Lorentz invariant but only Galilean 



invariant. 

The large N solution of the TDGL equation leads to the following asymptotic expressions 



for the effective mass and correlation function [30] 



M 2 (t) t =°° A (6.6) 
4r 

(rj a (z, r) V b (6, r)> = 6 a > b S exp (6.7) 

where again the averages are over the initial probability distribution and S Q depends on the 
initial conditions. 

Furthermore, recently a detailed study of the asymptotic behavior in time in the large iV 
limit of phase ordering kinetics in these models revealed a form of Bose Einstein condensation 



of long wavelength fluctuations asymptotically in time |3T . 

The TDGL description leads asymptotically to Goldstone bosons, with an effective mass 
term that vanishes asymptotically as in eq.(|6.6|), whereas in QFT the effective mass vanishes 



as in eq. ( 4.16 ). The oscillatory behavior of ( |4.16| ) is an important difference because it is 



responsible for non-perturbative anomalous dimensions |[L1J| . The dynamical correlation 
length obtained from the TDGL equal time correlation function, eq. ( |6.7| ) agrees (up to an 
overall scale) with that obtained in QFT for r < t s and given by eq.( |5.4| ). However, whereas 
the TDGL correlation length holds at all times (after some short initial transients [^0| ,^l|] ) , 
the QFT solution exhibits a crossover at r > r s [see eq. (|5.14 )1 to a linear dependence 



determined by causality and the onset of a Bose-Einstein condensate. This crucially different 
behavior is related to the non-relativistic nature of the TDGL equations in contrast with 
the fully relativistic invariance of our model. 

A remarkable similarity is that both the QFT version and the TDGL one of the large 
N model dsiplay the onset of a novel form of a Bose Einstein condensate. In QFT we have 
seen that the q = quantum mode grows linearly with time as a consequence of a vanishing 
effective mass term. This condensation mechanism is responsible for cutting off the spatial 
correlations for r > It as discussed above, a result compatible with causality. In the TDGL 
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description, a Bose condensate is found []3T| in order to satisfy the consistency condition in 
the asymptotic equilibrium state, this condensate results in a Bragg peak at zero spatial 
momentum ||31|| . An important difference is that whereas the TDGL description is purely 



dissipative and the asymptotic dynamics is rather insensitive to the initial conditions and 
coupling, in the QFT description dissipation is compatible with energy conservation and 
time reversal invariance and the asymptotic evolution depends on the initial conditions and 
coupling. 



B. Textures and scaling in a FRW background: 



In references [32| the process of phase ordering in radiation and matter dominated FRW 
cosmological backgrounds has been studied at the classical level in the non-linear sigma 
model (NLSM). A more detailed connection between the cosmological setting and the TDGL 
description described above has been provided in reference |3B| and we refer the reader to 
those references for more details. We here highlight the most important ingredients of the 
solution to compare to our results on quantum phase ordering kinetics which obviously are 
valid for Minkowski space-time. The main evolution equation studied in refs. [|3~2| , |3~3| are 
given by (using the same notation for the variables) 



d 2 fj ot(r) dfj 



dr 2 



dr 



V 2 ff+M 2 (r) rf= 



(6.8) 



where r now refers to conformal time, ot(r) = 2 d In a(r)/d In r and a(r) the FRW scale 
factor as a function of conformal time. Here the effective mass Ai 2 (r) arises from a spatial 
Gaussian average of the non-linear term in the equation of motion of the non-linear sigma 
model 



33|. The scaling solution found in references [q2|,q31] result in the self-consistent 



behavior for the effective mass term given by 

MHt) = 



(6.9) 



which leads to an exact form for the solution of the mode functions (in terms of comoving 



wavevectors) in terms of Bessel functions |32|,|33|]. The solutions are of the scaling form 



%(t) <x F(\k\r) 



(6.10) 



As shown in |J3] these gaussian scaling solutions lead to equal (conformal) time correla- 
tion functions that are cutoff at r = 2r and that reveal a dynamical (comoving) correlation 
length that grows linearly in conformal time £ ~ r. In radiation and matter dominated 
FRW cosmologies the size of the causal horizon ~ a(r) r in conformal time and the cutoff 
in the correlation function is a consequence of causality. 

The solution of refs. [ |3"2"| , |3"3"| have no a = limit. The normalization and the power 
spectrum have an infrared divergence for a — > which makes impossible to fulfil the NLSM 
constraint in that limit. Thus, our quantum Minkowski solution is not related to the classical 
gaussian solutions in FRW for the NLSM. Furthermore, our treatment is fully quantum 
whereas a classical stochastic approach is used in refs. p2|,|33 |. 
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There are only two aspects where our quantum results and the solution of refs. p2|j33 



are similar: a correlation length that grows linearly with time and a correlation function 
with support inside r < 2r. 

We find a crossover from the spinodal regime to the non-linear regime. Of course this 
crossover cannot be captured by the NLSM since in the NLSM the field is constrained to lie 
at the minima of the potential, i.e. there is no possibility of spinodal instabilities or rolling 
of the expectation value. In addition, our effective mass squared decreases as 1/r while it 
oscillates with frequency 2 m whereas in the NLSM the effective squared mass vanishes as 
1/r 2 . This different behaviour corresponds to the different nature of the IR divergences in 
FRM and Minkowski space-times. 



VII. CONCLUSIONS 

In this article we have obtained a novel set of non-perturbative and non-equilibrium phe- 
nomena in the leading order in the large N limit of a scalar theory. We have elucidated the 
shortcomings of a static effective potential description in the large N limit and contrasted 
this formulation to the results obtained within a consistent and non-perturbative treatment 
of the dynamics out of equilibrium in the same approximation. Whereas the static effective 
potential is complex in the region of coexistence, the non-equilibrium evolution leads to a 
dynamical Maxwell construction: all the expectation values in the coexistence region are 
allowed asymptotic solutions of the dynamics and the "restoring force" for the evolution of 
the expectation value vanishes asymptotically, i.e. a flat dynamical potential in the coex- 
istence region. The dynamical Maxwell construction is a consequence of non-perturbative 
particle production through spinodal instabilities. The particles are produced mainly at low 
momentum. A dynamical correlation length emerges whose behavior as a function of 
time reveals a crossover from oc \fi for t < t s ~ |m| _1 ln(l/A) to oc (£ — t s ) for 
t > t s . For t > t s the evolution is dominated by the non-linearities and a true scaling solu- 
tion emerges for which the field scales with a non-perturbatively large anomalous dynamical 
exponent z = 1/2. These novel phenomena in the non- linear regime are a consequence of 
the onset of a non-equilibrium Bose-Einstein condensate at zero momentum that grows lin- 
early in time. This condensate and the density of modes nearby ensure that the equal time 
two-point correlation function vanishes by causality for r > 2(t — t s ). We have studied the 
process of quantum phase ordering and interpreted the long-time phase ordering dynamics 
as proceeding via the growth of domains or bubbles that expand at the speed of light inside 
which there is a non-perturbative condensate of Goldstone bosons. The equal time two-point 
correlation function falls off as 1/r inside these domains. 

The analysis in terms of a quantum density matrix leads to a description of the non- 
equilibrium ensembles and correlations in terms of a classical but stochastic distribution 
function which describes large amplitude configurations with unsuppressed probabilities. 
We have also compared our results of quantum phase ordering dynamics to those obtained 
in stochastic field theories both in condensed matter (non-conserving TDGL equation) as 
well as the non-linear sigma model in a Friedmann- Robertson- Walker cosmology establishing 
some qualitative similarities and important differences. 

Furthemore we have also clarified aspects of the dynamics when the initial energy density 
is larger than that of the top of the potential. Our analysis shows how the expectation 
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value is driven to zero, showing unambiguously that the dynamics is O(N) symmetric as a 
consequence of non-perturbative particle production via parametric amplification. 
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Fig. 4: eta(tau) for g=10 A -5, eta(0)=5 (mu A 2 = -1 < 0) 
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